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Abstract. We study the effects of correlated low frequency noise sources acting 
on a two qubit gate in a fixed coupling scheme. A phcnomenological model for 
the spatial and cross-talk correlations is introduced. The decoherence inside the 
SWAP subspace is analyzed by combining analytic results based on the adiabatic 
approximation and numerical simulations. Results critically depend on amplitude of 
the low frequency noise with respect to the qubits coupling strength. Correlations 
between noise sources induce qualitative different behaviors depending on the values 
of the above parameters. The possibility to reduce dephasing due to correlated low 
frequency noise by a recalibration protocol is discussed. 
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1. Introduction 

Solid state nanodevices are at the forefront of present research towards the 
implementation of quantum networks for quantum computation and communication. 
The impressive development in device design and control tools achieved in recent years 
has by now to face intrinsic limitations due to material imperfections and fluctuations. 
The resulting noise presents a variety of material and device dependent features, 
ranging from noise spectra showing narrow resonances at selected frequencies (sometimes 
resonant with the nanodevice relevant energy scales) to low-frequency high-amplitude 
noise, often displaying 1/f behavior. Modelling these fluctuations has naturally lead 
to overtake the ubiquitous effective bath description via harmonic models and/or the 
hypothesis of linear coupling to the device under investigation. 

A typical example are background charge fluctuations which are known since more 
than 10 years to strongly affect the performance of Single-Electron Tunneling (SET) 
circuits pQ. Nowadays they represent the main limitation for any nanocircuit gate 
requiring highly reliable electrostatic control. This is clearly the case of charge [2] 
and charge-phase [3] superconducting qubits, but also of semiconducting spin qubits 
elecrostatically coupled to form a two-qubit gate pE]. 

The general belief is that background charge noise is due to the activity of random 
traps for single electrons in dielectric materials surrounding the island of SET devices or 
of superconducting nanocircuits. These traps may have different trapping energies and 
switching times, 7 _1 . An ensemble of non interacting traps with a uniform distribution of 
trapping energies and a I/7 distribution of switching rates may originate the frequently 
observed 1/f noise [5]. Such a spectrum is indicative of numerous traps participating 
in the generation of the noise. On the other hand, some samples clearly produce a 
telegraph noise with random switching between a few states (with a magnitude of up to 
0.1 e in SET devices) [HE]. In addition, recent observations on superconducting qubits 
in different setups, have suggested the possibility that a few impurities may entangle 
with the device [8] . Such a variety of experimental facts may be consistently predicted 
by describing background charges as two-state systems whose dynamical behavior may 
turn from quantum mechanical to classical with increasing temperature and/or with 
increasing the strength of their dissipative interaction with the fluctuations of the 
surrounding local host [9] . Such a modellization clearly departs from the " conventional" 
bosonic bath model and nicely predicts the different protocol-dependent decay laws 
of the coherent dynamics observed in charge and charge-phase qubits [9\. Multiple 
frequencies in the qubit dynamics and dependence on the uncontrollable impurities 
initial state at the beginning of measurement protocol are typical manifestation of non- 
gaussian character of background charge fluctuations [10J. 

In order to limit the effect on single qubit gates of these material-specific 
fluctuations different strategies have been developed. Amongst the most successful 
is the design of nanocircuits operating at "protected working points" insensitive to 
charge fluctuations to lowest order in the noise strength [3]. Open- and closed- loop 
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control protocols partly mutated from quantum optics and NMR [TT] , represent another 
promising route [12], [131 EH] ■ 

Presently, the effort of the scientific community working with Josephson qubits is to 
appropriately extend the above strategies to multiqubit architectures, the first step being 
to implement an efficient two-qubit gate. Different schemes to couple superconducting 
qubits have been proposed [15] and some experiments pointed out the possibility to 
realize the desired entangled dynamics [T6| [T71 [T8 J. However, achievement of the 
needed high fidelity is still an ambitious task. In addition to fluctuations experienced 
individually by each single qubit gate, coupled qubits, being usually built on-chip, may 
suffer from correlated noise due to sources acting simultaneously on both sub-units. 
The effects on two-qubit gates of uncorrelated and correlated bosonic baths has been 
investigated [T9] . 

As far as the effect of background charges, opinions about the probable location 
of traps are divided and observations depend to a certain extent on the specific sample 
and on the junction geometry. However, there is unambiguous evidence that fluctuating 
traps located in the insulating substrate contribute essentially to the total noise in SET 
devices [UEO]- These traps are expected to induce similar fluctuations on the two islands 
built on the same substrate. On the other hand fluctuating traps concentrated inside 
the oxide layer of the tunnel junctions, due to screening by the junction electrodes, are 
expected to act independently on the two qubits [I]. 

Fluctuating impurities acting simultaneously on coupled qubits represent a further 
unconventional noise source which solid state nanodevices has to face [21] . This is the 
subject of the present paper. Specifically, we will introduce a model for correlated charge 
noise on interacting charge-phase qubits in a fixed capacitive coupling scheme. Relying 
on measurements on SET circuits of power spectra on the two transistors and of the 
cross-spectrum power density, we suppose a 1/f behavior for both the two channels 
spectra and the cross-spectrum pQ. In addition the cross-talk between the two qubits 
due to the capacitive coupling itself between the islands will be discussed. Our analysis 
is based on analytical results obtained within the adiabatic approximation for the 1/f 
noise and on the numerical solution of the stochastic Schrodinger equation. Solving the 
dynamics from short-to-intermediate time scales allows complete understanding of the 
effects of correlations. Our work extends the analysis of Ref. [2 1 J which, being limited to 
the long-times behavior, misses relevant features occurring in the short time domain. We 
find that usually correlations induce a faster decay of the coherent dynamics compared 
to the action of independent fluctuations. Nevertheless, under realistic values of low- 
frequency noise amplitude, increasing the degree of correlation may instead lead to 
longer decoherence times. Finally, the possibility to reduce the effects of low-frequency 
correlated noise via open-loop recalibration protocols is discussed. 

The paper is organized as follows: in Section 2 we introduce the setup consisting 
in two Cooper pair boxes coupled by a capacitor. The cross-talk effect and the charge 
noise sources responsible for correlations will be described and their correlation factor 
defined. In Section 3 we present a possible model for correlated noise exhibiting 1/f 
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power spectrum and cross-spectrum. In Section 4 relevant dynamical quantities for a 
two-qubit gate are introduced and analytic/numerical methods are illustrated. Section 
5 includes our results for the entangled qubits dynamics in the presence of correlations. 
Conclusions are drawn in Section 6. 



2. Coupled Cooper-pair boxes and noise correlations 

In the fixed capacitive coupling scheme for charge [2] or charge-phase |3j qubits the 
islands of two Cooper pair boxes (CPB) are connected through a capacitance [HH E]> 
as illustrated in figure HJ The system is described by the Hamiltonian 

7-Lo= ^2 ^a + ^CcO?l-?l,x)(<72-?2,x)> i 1 ) 
ae{l,2} 

where each CPB is modelled by 

T-L a = [E a! c{q a - q a , x ) 2 + E a j cos0 a ] . (2) 

E a ,c = 2e 2 /C Q , i s is the charging energy of the island belonging to CPB a, the 
total island capacitance C a ^ = C a ^G + C a ,j being the sum of the gate and junction 
capacitances. q a>x = C a ^V a ^ / {2e) is the corresponding dimensionless gate charge. 
Cooper pair tunneling across the Josephson junction a requires an energy E a j. E cc = 
(2e) 2 C T /(C , i iS C , 2 ,s) is the coupling energy, with 1/C T = 1/C C + + 1 / C 2,y, the 

total inverse capacitance of the device. 
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Figure 1. Two CPB connected by a fixed capacitor. 

The dimensionless charge q and the phase (p of each box are conjugated variables, 
[<p, q] = i. The system is subject to fluctuations of different origin. In part they arise 
from the control circuitry and can be described by an effective impedance modelled by a 
conventional bosonic bath. Noise sources of microscopic origin are atomic defects located 
in the oxide of the tunnel junctions, leading to fluctuations of the Josephson energy and 
background charges acting like additional uncontrollable q a ^ sources. Devices based on 
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the charge variable are particularly sensitive to background charge fluctuations. Usually 
they can be modelled as two-level fluctuators (TLF) inducing a bistable polarization of 
the superconducting island. A collection of TLF produces a noise whose spectral density 
approximately follows a 1/f law. They are responsible for the sensitive initial reduction 
of the amplitude of coherent oscillations in single qubit gates observed when repeated 
measurements are performed [9j [12] . Fluctuations of polarization islands are expected 
to be a severe hindrance for coupled qubits gates based on the charge variable [22] . Low 
frequency charge fluctuations lead to an additional stray contribution to the gate charge 
g QjX , which can be modelled by a random variable Sq aiX (t) leading to 



Note that the coupling capacitance induces a cross-talk between the two devices, i.e. 
fluctuations <5gi )X (t) acting on q 2 and vice-versa. As we already mentioned, background 
charges responsible for gate charge fluctuations are spatially distributed in a device- 
dependent unpredictable way. Possibly they are partly located in the substrate, partly 
in the oxide layer covering all electrodes, partly in the oxide barriers of the tunnel 
junctions. Due to the shielding by the electrodes, impurities within tunnel junction a 
are expected to induce only gate charge fluctuations 5q a , x (t). On the other hand, random 
arrangement of noise sources in the bulk substrate originate correlations between gate 
charge fluctuations to an extent depending on their precise location [TJ. Pictorially, 



Figure 2. Pictorial representation of cross-talk (a) and spatial correlations (b). (a) 
The uncorrelated noise sources s\ and S2 act on coupled qubits: because of the coupling 
each qubit also suffers from the noise source directly acting on the other qubit. (b) 
Two non interacting qubits in the presence of s\, S2 and of S12 that simultaneously 
acts on both qubits. s%2 may represent a set of of impurities located in the insulating 
substrate. It generates a fluctuating interaction even in the absence of direct coupling 
between the two qubits. 



we may separate impurities in two ensembles {si,S2} and s±2 influencing either each 
sub-unit or both islands, as illustrated in figure [2j The noise 5q a ^ viewed by qubit a is 
due both to s a and S\2- Correlations between 5q\ yX and Sq2, x originate from set s\2, are 
termed spatial correlations. 

In the above phenomenological description, we assume that $q a ,x{t) are stationary 
stochastic processes having zero average and the same variance a 2 . We quantify the 
degree of spatial correlation between £g l x and 5g 2 ,x via (&Za, x (£)&2 , /3,x(£)) = [<W + /%>(! — 
8 a 0)]a 2 . Evaluation of /i sp would require a microscopic description of the device and it 



7i = Ho + 57i, 
5H = 



(3) 



qi[2E ltC Sqx, x (t) + E cc 5q^ x (t)} - q 2 [ E ccSqi, x (t) + 2£ 2 , c 5g 2 ,x(*)] 
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is expected to depend on the dimension of the boxes and on their relative distance, on 
the specific spatial distribution of impurities s 12 and on the relative weights on 5q a ^{t) 
of fluctuations due to set s\ 2 and s a , as shown in Ref. pQ. 

Cross-talk due to coupling and spatial correlations enter the overall noise X a felt 
by each sub-unit: 

Xt{t) = 2E 1>C 5q 1<x (t) + E cc Sq 2>x {t) (4) 

X 2 {t) = E cc Sqi, x {t) + 2E 2 , C 6q 2> *{t) (5) 

The amount of correlation between X\ and X 2 may be quantified by a correlation 
coefficient, which is defined for general stochastic processes £i and £2 as [24"] : 



= (mi - yfeg) - jj) (6) 

where (•) indicates the ensemble average, and ^ = (&(£)). ^From (J4]) and (jSJ) we obtain 
(X*> = (4i% c + e 2 cc + 4/i sp £ a , c £ cc )a 2 , 
(X a X 2 ) = [2E CC {E 1>C + #2,0) + /is P (4£i, c £2,c + ^cc)]^ 2 

thus the correlation coefficient of X\ and X 2 reads 

^ 2-E cc (-Ei ;C + E 2yC ) + fJ> sp (4:E^cE 2) c + E cc ) 

^El c + E 2 CC + 4/i sp £ 1 , c £cc)(4£ 2 2 c + E 2 CC + 4 yUsp £ 2 , c £ C c) 

In the absence of spatial correlations 5q a ,x are independent and only the effect of cross- 
talk is left. The correlation coefficient in this case reduces to 

^ ^EccjE^c + E 2 ,c) ^ 4E c E C c ^ 

'{AE 2 + E GC )(AE 2 + E cc ) ~ ( 4 ^c + E co ) ■ 



In the last approximation we have supposed E a c — Eq, within manufacture tolerances. 
For typical values of parameters for charge qubits E c ^> E C c thus /i ct ~ E C c/Ec, 
giving values between 0.015 [22J and 0.12 [T7]. Clearly, larger values of the coupling 
strength E C c, desirable to produce faster two-qubit gates, would also lead to higher 
cross-talk correlations /i ct . In general, the correlation coefficient (EJ) for E a> c — Ec is 
approximately given by 

''' ~ AE 2 C + E 2 CC + Afx sp E c Ecc ~ 1 + /Ws P 
where we used OH]). Strong correlations between Xi and X 2 , fim 1, may originate either 
from large cross-talk or from large spatial correlations. For instance by engeneering 
device design [23] it could be possible to get /i ct — 1 implying /i ~ 1. On the other 
hand, in the presence of a low level of correlations /x sp /x c t ^ 1, equation (|U]) simplifies 

tO yU ~ ii sp + /i ct . 

In general ([7]) gives the overall amount of correlation between fluctuations affecting 
the two CPBs. In the following we will not specify the physical mechanism responsible 
for these correlations. We will simply suppose the existence of a degree of correlation 
between the fluctuations X\ and X 2 quantified by the coefficient /x. 
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3. Charge noise power spectra and cross-spectrum 

Measurements of charge noise due to background charge fluctuations in SET devices pQ 
have revealed a 1/f behavior at low frequencies (measurements extend down to about 1 
Hz), with a roll-off frequency of 100 — 1000 Hz. Two SET whose islands are positioned 
about 100 nm apart show a similar 1/f behavior for the cross-spectrum (defined later 
by (TTTT) ) indicating correlations between fluctuations affecting both islands. Similarly, 
measurements of charge noise in charge-phase [12] qubits shows a 1/f behavior for 
/ < 100 kHz whose amplitude depends on temperature, on junction size and on screening 
of the island by electrodes. Echo experiments suggest that 1/f noise extends up to 1 
MHz. In this setup charge noise at higher frequencies (up to 10 MHz) is due to driving 
and readout subcircuits and it is characterized by white spectrum. Measurements of 
energy relaxation processes in charge qubits have suggested that charged impurities 
may also be responsible for ohmic noise at GHz frequencies [25]. To our knowledge, 
measurements of cross-spectrum on these class of nanodevices have not been reported 
in the literature. It is however expected that, similarly to SET devices, correlations 
between fluctuations acting on superconducting islands of the two on-chip CPBs display 
1/f cross-spectrum at low frequencies. 

Our goal is to introduce a model for the fluctuations X\(f) and X 2 (t) such that 
both power spectra and cross-spectrum 

S Xa (u) = / °°dre- w [(X a (t + r)X Q (t)) -X 2 a ] (10) 

J — oo 

/+oo 
dre-^ T [(X 1 (t + r)X 2 (t))-X 1 X 2 ] (11) 
-oo 

display similar 1/f behavior at low frequencies and are characterized by a finite 
correlation coefficient defined by (jSJ). To this end we introduce two independent 
stationary stochastic processes, n\(t) and ri2(t) with the same average and characterized 
by the same autocovariance function and spectrum [23] 

C na n a (r) = (n a (t + T)n a (t)) -n\ = C(r) 

r+oo 

Sn a nM= / dTe-^ T C(r)= S(u). (12) 



The processes Xi(t) and X 2 (t), defined as linear combinations of nx(t) and n 2 {t) 
X 1 (t) = y/l- V n 1 (t) + Vvn 2 (t) 

X 2 (t) = Vvni(t) + ^/l~^n 2 {t) , (13) 

with i] E [0, |], are correlated and their correlation coefficient reads 

/x = 2^/77(1 - 77) . (14) 

Thus fj, is a monotonic function of r\ £ [0, |] ranging in the interval [0, 1]. If rj = 0, X\ 
and X 2 reduce respectively to the uncorrelated processes ri\ and n 2 and \i = 0. Instead 
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when V — k the correlation factor reaches its maximum value \x = 1. In this case X\ 
and X 2 reduce to the same process: 

/i = 1 => X^t) = X 2 (t) = ^=[ni(t) + n 2 (t)] . (15) 

The autocovariance functions of Xi(t) and X 2 (t) are identical and read 

C Xa x a (r) = (X Q (t+T)X Q (t))-Xl = (I-77) Cn ini +vC n2 n 2 = C(r) , (16) 

therefore Xi and X 2 have the same variance a 2 and equal power spectra S Xi {uj) = 
Sxi{w) = S(u) given by (jT2J). This simple model for correlated noise allows, by 
changing the arbitrary parameter rj, to modulate the correlation coefficient between 
X\ and X 2 from to 1, maintaining the desired spectrum S{u) for both processes. 
It is worth noticing that the first order statistics of X% and X 2 depends on 77: 
X\ = X 2 = n{\/l — rj + ^yfj). To avoid this dependence we set n — 0, implying vanishing 
average values for X\ and X 2 . The correlation factor enters the cross- covariance and 
the cross-spectrum of X\ and X 2 [21] 

C XlX2 (r) = {X 1 (t + r)X 2 (t))-X l X 2 = hC(t), (17) 

/+00 
dre-^ T C XlX2 (r) = fiS(u). (18) 
-00 

It can therefore be detected by spectral analysis via [26] 

=At . (19) 

y/S Xl (uj)S Xi (uj) 

By measuring power spectra and cross-spectrum of voltage fluctuations across each SET 
in the frequency range 1 to 10 Hz, Zorin et al. estimated according to f|T9|) the correlation 
coefficient /1 = 0.15 ± 0.05 pp. 

In order to obtain a 1/f spectrum for processes X a (t) we adopt a commonly 
employed model which consists of an ensemble of independent TLF. Each fluctuator 
incoherently switches between two metastable levels, with a rate 7^, producing a random 
signal £fc(i). This signal has a lorentzian power spectrum, S^ k (uj) = hvljk/(lk + Co ' 2 )' 
Vk being the difference between the values assumed by £k{t)- When the switching rates 
jk are distributed according to P(j) oc I/7 in [7 m , 7m], the overall noise obtained by 
summing all TLFs contributions displays al/w behavior in [7 m , 7m] [H] 

{ =E&w - ^) = E^^)^ (20) 

/c fc=i v y 

where A = ti(v 2 )Ntlf/[^ ln(7M/7m)] and A^tlf is the total number of fluctuators. If 
the independent random processes n\(t) and n 2 (t) are generated as a sum of such an 
ensemble of TLFs, the spectrum of each n a (t) will be 1/f in [7 m , 7m] and will have 
variance a 2 = 5- / duS(u) = \Ntlf(v 2 ) ■ The 1/f correlated stochastic processes X\(t) 
and X 2 {t) are obtained from f|T3|) once the phenomenological correlation factor /1 is fixed. 
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4. Two qubit gate and relevant dynamical quantities 

At sufficiently low temperatures each CPB may operate as an effective two-state system, 
the coupled boxes implementing a two qubit gate. In the fixed coupling scheme the 
interaction is switched on by individually manipulating each qubit to enforce mutual 
resonance conditions p2J [17]. This allows the realization of elementary two qubit 
operations. We denote the lowest eigenstates of each CPB as {|+) Q , |— ) a }, with splitting 
depending on the control parameter q a>x , Q a (q a ,x)- By operating at the so called "charge 
protected point", q a>x = 1/2, the system is insensitive to charge fluctuations at lowest 
order, meaning that dQ a (q a}X ) / dq a , x \q a ^=o.5 = [3J . In a pseudo-spin description, 
in the eigenstate basis charge fluctuations are off-diagonal at this working point. 
Therefore projecting the coupled boxes Hamiltonian (j2J) into the computational subspace 
® Ij'M - i,j e {+, -} we get 

H = Hq + 5H (21) 

Ho = ® n (2) - ® 4 2) + ® -S 2) (22) 

5H = -^®1^ - ^®af\ (23) 

where we assume the two qubits are tuned at the same Bohr splitting Q, whose typical 
value is of the order 10 11 rad/s. Here E C c = 2£?cc Qi,+- ?2,+- an d X<* — 2A Q g Q , + _, 
being = (+\ a q a \-)a- 

We remark that this symmetric configuration is hardly reachable in practice by 
fabrication accuracy only. In the charge-phase implementation [3] this can be achieved 
thanks to the characteristic two-port design. The quantization is in fact based on a 
split CPB connected to a large measurement Josephson junction, the phase across it, 5, 
representing an additional control knob. The device presents a "doubly" protected point 
at q x = 1/2 and 5 = 0, which is a saddle point of the single-qubit energy splitting versus 
external parameters. In the two-qubit setup resonance is achieved by slightly displacing 
one of the qubits from the phase protected point, but maintaining both qubits at the 
charge protected point. This setup is therefore expected to be sensitive also to phase 
fluctuations and this topic will be addressed elsewhere [27] . 

The setup described by (1211) may in principle implement a iSWAP gate. This is 
easily illustrated in its eigenstate basis, reported in tabled] in terms of the dimensionless 
coupling strength g = Eqc/VL. In the absence of fluctuations, the two-qubit Hilbert 
space is factorized in two subspaces spanned by pairs of computational states. In 
particular, the system prepared in the state \-\ — ), freely evolves inside the subspace 
spanned by {\-\ — ), | — h)}, reaching the entangled state (|H — ) + i| — h))/v2 at time 
f = iQ = n/2g. States {|1), |2)} generate the so called SWAP subspace, whereas we 
refer to Z subspace as the one spanned by {|0), |3)} [27] . 

A numerical analysis has shown that for typical values of parameters in charge- 
phase qubits, the SWAP-eigenvalues are more stable than single qubit splitting with 
respect to uncorrelated gate charge fluctuations [28]. As a consequence, for sufficiently 
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Table 1. Eigenvalues and eigenvectors of Hq. Here sin^ = g/ (2a/ 1 + g 2 /4) and 
cos<^ = + g 2 /4. 



i 


A (0 


I*) 





-Vi + .g 2 /4 


[-(sin V /2)|++) + (cos V /2)|--)] 


1 


~5/2 


H + -> + |-+»/v^ 


2 


5/2 


(|+-) + |-+))/V2 


3 


V1 + .9V4 


[(cos^/2)| + +) + (sin^/2)|--)] 



small-amplitude low-frequency charge fluctuations the decay time for entangled states in 
the SWAP subspace is expected to be longer than the single qubit dephasing time [28J. 
In the following we analyze how low frequency correlations between charge noise felt by 
each qubit may influence the dynamics in the SWAP subspace. 

4-1- Relevant dynamical quantities 

In the presence of low frequency fluctuations the calibration of the device is unstable. 
As a result, the quantum dynamics of the interacting qubits will depend on the 
measurement protocol, as already observed for single qubit gates. Ideal quantum 
protocols assume measurements of individual members of an ensemble of identical 
(meaning that the preparation is controlled) evolutions, defocusing occurring during 
time evolution. In practice for solid state nanodevices several samples are collected 
during an overall measurement time t m . Lack of control of the environment preparation 
determines defocusing of the signal, analogous to inhomogeneous broadening in NMR. 
The considerable initial reduction of the amplitude of the coherent oscillations of single 
qubit gates affected by 1/f charge noise is due precisely to this effect [9j [12]. On 
the other hand the effect of low frequency noise on relaxation processes is negligible. 
Thus the system dynamics can be treated in the adiabatic approximation for the low 
frequency charge fluctuations. Under this approximation scheme populations of the 
system eigenstates do not evolve. The relevant dynamical quantities are therefore the 
off-diagonal elements of the system density matrix in the same basis. 

The efficiency of the iSWAP protocol in the presence of 1/f spectra on both qubits 
and 1/f cross-spectrum can be therefore extracted by evaluating a single dynamical 
quantity, the coherence between the eigenstates of the SWAP subspace. The two 
qubit density matrix in the presence of the dimensionless classical stochastic processes 
Xi(t) = Xi(t)/Q generally reads 

p(r) = J P[a; 1 (r')]P[x 2 (r')]P[a: 1 (r'),a; 2 (r')]p[r|a: 1 (r'),x 2 (r')], (24) 

where p[r|xi(r'), x 2 (t')] is the system density matrix calculated for a given realization 
{xi(t'), x 2 (r')}. The integration is over over all possible realizations weighted by the 
probability density P[xi(r), x 2 (r)]. 
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We solved equation (j24p numerically by generating the independent random 
processes U\{t) and n 2 (r) and from them the correlated processes Xi(t) and x 2 (t), as 
illustrated in section [3j The Schrodinger equation related to the Hamiltonian ( 12T]) is 
numerically solved by a fourth order Runge-Kutta algorithm |29j, calculating the system 
dynamics p[r\xi(r), x 2 {r)]. These operations are repeated to perform an average of 
p[t\xi(t), x 2 (t)\ over many (> 10 4 ) realizations of the stochastic processes. Numerical 
simulations confirm that in the presence of low frequency noise (with 7m < 10 -2 Q) 
transitions between the SWAP and Z subspaces can be neglected. This further 
legitimates focusing on the coherence in the SWAP subspace which reads 

(l|p(r)|2> = Pl2 (r) = p 12 (0)e-^e- i$ M = p 12 (0) 

• J V[x 1 (r')}V[x 2 (r , )}P[x 1 (r'),x 2 (T , )]exp 

where u>i 2 [xi(t"),x 2 (t")] gives the noise renormalized splitting between states |1) and 
1 2). The imaginary part of $(t) describes the decay of the entangled dynamics in the 
presence of adiabatic correlated noise. Further insight can be obtained by approximating 
(125]) to include the dominant inhomogeneous broadening effect. This is performed by 
applying the static path approximation (SPA), x a (t) = x a , which accounts for the lack 
of control of the device calibration via a statistical distributed gate charge at each run 
of the measurement protocol. In the SPA the coherence (125|) reduces to the evaluation 
of an ordinary two-variables integral 

Pu{r) =pvi(0) J dx 1 dx 2 P(x 1 ,x 2 ) exp[iruj 12 (x 1 ,x 2 )) , (26) 

where P(xi,x 2 ) is the joint probability density function of the random variables 
x\ and x 2 [24] . In the following we will use the notation (f(xi,x 2 )) to indicate 
/ dxi dx 2 P{x u x 2 ) f(x u x 2 ). 

In the following Section we will analytically evaluate the coherence in the SWAP 
subspace within the SPA in selected parameter regimes where numerical simulations 
have confirmed its accuracy. A numerical analysis will be performed to estimate the 
decay of entanglement under more general conditions. 



i / <1t"u 12 [x 1 (t" 



),X 2 (T 



")] 



,(25) 



5. Dephasing in the SWAP subspace: effects of correlations 



The average in (1261) is conveniently evaluated by performing the change of variables (fl~3]) . 
In fact, since the independent random processes n\(t) and n 2 (t) are generated from a 
large number of TLFs, their initial values n a are Gaussian distributed 

Pin.) = -i- exp[-n 2 /(2a 2 )] . (27) 



2trx 

Clearly, x\ and x 2 are two correlated Gaussian variables whose joint probability density 
function is (for \p\ < 1) [24] : 



P(x 1 ,x 2 ) 



1 



: exp 



1 



2a 2 (l - /i 2 ) 



[x\ + x\ — 2px 1 x 2 ) 



(28) 
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Figure 3. Dephasing factor given by (|31| for different values of the correlation 
coefficient \i for "weak" amplitude noise, a = 0.02 < g = 0.1. The dashed grey 
curve refers to the uncorrelated case /x = 0, the dashed black to /x = 0.50, the thick 
black refers to /x = 1.00. The crossing of each curve with the dotted horizontal line 
at e _1 identifies to the estimated dephasing time relative to each value of /x, 72 (/x). 
Intersections with the curve corresponding to uncorrelated noise (/x = 0) identifies 
the time t*(u) < t%. Inset: enlargement for short times. The thick grey line 
corresponds to /x = 0.75. The validity of the SPA approximation has been checked 
against numerical simulations for stochastic processes exhibiting a 1/f power spectrum 
in a range [7 m ,7\i] = [1, 10 6 ] s _1 (not shown). 



The effective splitting in the SWAP subspace in the presence of charge fluctuations 
entering the average (J26l) may be evaluated by exact diagonalization of the Hamiltonian 
( I2T1) . The solution of the resulting fourth order polynomial is rather lengthy so we do 
not report it here. Relevant features can be extracted by expanding the splitting up to 
fourth order in x a and keeping the dominant terms in the coupling strength g j27j 

^12(^1, x 2 ) w g - ^-{x\ + x\) + — {x\ - x\) 2 = g + 5ui 2 (xi, x 2 ) . (29) 
A tig 

This expansion suggests that the system behavior depends on the relative weight of the 
amplitude of the noise, measured in the SPA by the standard deviation a entering ([28]) . 
and the strength of the interaction between the qubits, g. In the following we consider 
separately the two regimes of "weak" and "strong" amplitude noise, a < g and a > g 
respectively. We remark that the nonmonotonous dependence of the splitting on x a 
may lead to almost degeneracy between the renormalized levels of the SWAP subspace. 
This effect may be relevant when the interplay of low and high frequency components 
is considered [27] . 

Weak amplitude noise a <^ g: In this regime ( 1291) can be approximated by keeping 
terms up to second order in a so that Scu^Zi, ^2) = — f ( x i + The average 

p 12 ( T ) =p 12 (0)e^(exp(-ir^r 2 )) (30) 

can be easily evaluated and leads, for the dephasing factor 

= [1 + (ga\l - ^)rf] " 1/4 x [l + {ga\l + ^)rf] ~ 1/4 . (31) 



Pl2{T) 
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Figure 4. Dephasing factor D^ 2 (r) for strong amplitude noise, g — 0.01 < a = 0.08. 
Curves correspond to different values of the correlation coefficient, from bottom to 
top fi £ {0,0.50,0.75,0.8,0.85,0.90,0.95,1.00}. Correlations improve the system 
performance despite the large noise amplitude. 



The dephasing factor factorizes into two contributions having the form of the decay of 
the single qubit coherence at protected point in the SPA [9] with standard deviations 
cr-y/l ± jj,. An analogous result has been found in [21]. Equation ( 13TT) is shown in 
figure [3] for different values of the correlation coefficient \i. For comparison the curve 
corresponding to independent noise sources acting on the two qubits is also reported, 
^ ) i2~°( r ) = [1 + {9°' 2t ) 2 } ■ Interestingly, at short times increasing the correlation 
coefficient induces a faster reduction of the amplitude of coherent oscillations in the 
SWAP subspace. This behavior crosses over to a regime where instead increasing 
the correlation coefficient slows down dephasing. This occurs at times larger than t* 
identified by the condition D^ 2 (t*) = D° 2 (r*), which gives r*(/i) = l/(ga 2 ^/l - /i 2 /2). 
The crossover takes place at times shorter than the dephasing time where D^fo) = e _1 

tM = 1 m « > r*(fi) . 32 

ga 2 {l -/i 2 ) 

We remark that, for quantum computing purposes it is crucial understanding the 
behavior at times shorter than the dephasing time. For instance fault-tolerant quantum 
computation (30j [31], i.e. implementing reliable quantum operations even in presence 
of errors, requires errors to be maintained below a small threshold (typically eth ~ 
10 _4 -rl0~ 6 ). The error of the iSWAP gate under investigation may be simply estimated 
as (in the adiabatic approximation) 

e = 1-Mp(tM = l-^ 2 (r) (33) 

being the iSWAP target state. This leads, at the dephasing time r 2 , to an error 
about 0.8. The correlation coefficient dependence of the dephasing time r 2 [2TJ does not 
reveal relevant features of the iSWAP gate operation occurring at initial time scales. 

Strong amplitude noise a ^> g: For larger amplitudes of the noise, the dephasing 
factor displays a completely different behavior. Numerical results are illustrated in 
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Figure 5. (a) Variance as a function /j,, for small amplitude noise a < g. (b) 
Correlation coefficient where the variance attains its minimum, as function of the ratio 
cr/g. For a > 4g it is fx m ~ 1. (c) Equation ([55)1 as a function /x for <r 3> g. 
decreases by different orders of magnitudes when /j, > 0.8. In (a) and (c) g = 0.01. 



figure HI In this case increasing the correlation coefficient systematically reduces 
dephasing. In particular the limiting case \x = 1, where the same noise simultaneously 
acts on both qubits, turns out to be relatively weakly affected by the low frequency 
noise, despite of its large amplitude. 

The short time behavior and its dependence on the correlation coefficient for weak 
and strong amplitude noise may be explained considering the renormalized splitting 
dependence on the noise variables x a . The level splitting itself in fact can be viewed 
as a random variable, with standard deviation as wi2 = \j {5uj\ 2 ) — {SoJ^p. The short 
times gaussian approximation for the dephasing factor gives 

DUt) * 1 - ^L 12 r 2 , (34) 

thus, at times short enough that a^ WV2 r < 1 (for data in figure El for r < 10 4 ), a larger 
deviation as^ 12 induces a larger dephasing, i.e. a smaller value for D% 2 (t). The variance 
of 5uu reads (see appendix for details) 

< 2 = ^{i(9 2 - ° 2 Y + A + A(9 2 + ^) 2 - 5a 4 ] + 2A 4 ]} • (35) 

For a < g, it reduces approximately to crj w „ = g 2 cr 4 (l + fi 2 ), monotonically increasing 
with fi, as shown in figure [5k. This explains the stronger dephasing observed for small 
amplitude noise at short times. On the other hand, for a > g, cr| £J12 is nonmonotonic 
with a minimum at /i m = [1 — |(cr/o) -2 — ^(cr/g) -4 ] 1 / 2 , which rapidly approaches 1 
(figure [5b). Thus for [l — > 1 the splitting variance rapidly reaches its minimum value, 
implying small dephasing at short times even for large noise amplitudes, a > 4g. 

5.1. Effects of higher frequencies and recalibration protocol 

We now consider the effect of correlated 1/f noise extending to higher frequencies 
maintaining the adiabaticity condition with respect to the qubit splitting Q, i.e. we 
simulate fluctuations leading to 1/f spectrum up to the cut-off frequency 7m = 10 9 s -1 . 
The resulting dephasing in this case, in addition to the inhomogeneous broadening 
mechanism, also originates from the dynamics of the fluctuators during the time 
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Figure 6. Effect of 1/f correlated noise extending in [7 m , 7m] = [10°, 10 9 ] s _1 for weak 
noise amplitude a < g = 0.1 and different values of fi. The dashed grey curve refers 
to [i = 0, the dashed black to [i = 0.50, the thick black to fi = 1.00. (a) Dephasing 
factor: correlations slightly increase dephasing at initial times and decrease it at longer 
times. The gray horizontal line refers to e _1 . (b) Effect of a recalibration protocol: 
increasing the correlation coefficient increases dephasing. Thin black line corresponds 
to fi = 0.25, thick gray to = 0.75. 



evolution. The dependence of the dephasing factor on the correlation coefficient jj, 
has the same characteristics observed in the presence of low frequency components only. 
The system still displays different behaviors depending on a being smaller or larger than 
g. This is illustrated in figure Eh and figure [TJi. In particular the very low dephasing 
already observed when the qubits are affected by the same environment (// = 1), persists 
also in the presence of high frequency noise components (figure UK)- 

In single qubit gates the inhomogeneous broadening effect may be sensibly reduced 
by a recalibration protocol resetting the initial value of the system polarization at each 
run of the measurement protocol [9]. Recalibration turns out to be effective on two- 
qubit gates also in the presence of correlations among the noise sources. Results shown 
in figure [6b and figure [7b have been obtained numerically by resetting the values of x a (0) 
at each realization of the stochastic processes x a (t). Interestingly, even if the effect of 
low frequency components is practically eliminated by the recalibration procedure, the 
decay has a different dependence on /i depending on a. In particular, if a < g the larger 
is the correlation coefficient, the faster the signal decays (figure [6b), if instead o > g 
stronger correlations correspond to slower decay (figure [7b) . 

6. Conclusions 

In the present paper we have introduced a phenomenological model for 1/f correlated 
noise affecting a two-qubit gate in a fixed coupling scheme. Our analysis is based on 
analytical results obtained within the adiabatic approximation and on the numerical 
solution of the stochastic Schrodinger equation. 

Due to the nonmonotonicity of the renormalized splitting in the SWAP subspace, 
the entangled dynamics sensitively depends on the ratio ojg between the amplitude 
of the low frequency noise and the qubits coupling strength. For small amplitude 
noise, correlations increase dephasing at the relevant short times scales (smaller than 
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Figure 7. Effect of f/f correlated noise extending in [7 m ,7M] = [10°, 10 9 ] s -1 for 
strong noise amplitude a > g = 0.01 and different values of fi. (a) Dephasing factor 
for fi = 0: black curves (light black for slow noise, tick black for low and high frequency 
noise); fi = 1: gray (slow noise) and red (slow plus high frequency noise), (b) The 
effect of a recalibration protocol is shown by the dashed curves, [i = black, // = 0.75 
gray, fi = 1 red. In both protocols (a) and (b) increasing the correlation coefficient 
decreases dephasing. 



the dephasing time). On the other hand, under strong amplitude noise, an increasing 
degree of correlations between noise sources acting on the two qubits always leads to 
reduced dephasing. Our numerical analysis has shown that the above features hold true 
for adiabatic 1/f noise extending up to frequencies 10 9 s -1 about two order of magnitudes 
smaller that the qubit Bohr frequencies. 

We remark that the above results on the reduced dynamics in the SWAP subspace 
apply also to diferent two-qubit gates involving the same states, at least as long as 
the adiabatic approximation holds true. The performance of two qubit gates involving 
states of the Z subspace, like the c-NOT gate, might be reduced in view of the larger 
sensitivity of the Z subspace to low- frequency charge noise [28] . 

We have analyzed the possibility to reduce the effects of low-frequency correlated 
noise by a open-loop recalibration protocol of the two-qubit gate. Despite counteracting 
the inhomogeneous broadening effect, the efficiency of the protocol still depends on the 
value of o~/g, the maximum efficiency occurring for small amplitude uncorrelated noise 
(o~ < g and \i = 0), or for strong amplitude correlated noise (a > g and \i = 1). 

The observed reduced sensitivity of the SWAP subspace to correlated strong 
amplitude noise might suggest exploiting this subspace to reliably encode a single qubit, 
in the same spirit of Decoherence Free Subspaces (DFS) [32]. The possibility of avoiding 
errors due to correlated noise by encoding in DFS has indeed been recently discussed 
for superconducting qubits in [23}, [331 El] • This strategy rigorously applies to the pure 
dephasing regime where a DFS subspace exits for collective noise. In the situation 
analysed in the present article however the SWAP subspace is not rigorously a DFS. In 
fact, in the presence of collective noise (fj, = 1) the interaction Hamiltonian (1231) reads 
— ^(a^ (g> 1^ + (g> o~[ )X. The system operator entering this coupling term has two 
degenerate eigenstates which do not span the SWAP subspace. In addition, for finite 
values of g, the system Hamiltonian (1221) does not leave invariant the subspace spanned 
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by the degenerate eigenstates, as DFS should require [35J. 
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Appendix A. Moments entering o^ 12 

From ( |28|) it comes out that the marginal probability density function of x a is a Gaussian 
function with standard deviation a and zero average value. Then: 

{xln) = f^T 2n and ^ n+1 > = ' ^ 

see Ref. [36]. Evaluation of the variance of the splitting fluctuations as^ 12 requires 
knowledge of the following mixed moments: 

(xlxl) = ^(1 + 2/i 2 ) 
(x\x\) = 3a 6 (l + V) 
(xlx b 2 ) = 15a 8 (l + 6/i 2 ) 

{x\x\) = a 8 (9 + 72/i 2 + 24//) (A.2) 
which directly follow from (jlBl by using flA.ll) . 
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